The Little Rip 
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We examine models in which the dark energy density increases with time (so that the equation- 
of-state parameter w satisfies w < —1), but w — s> — 1 asymptotically, such that there is no future 
singularity. We refine previous calculations to determine the conditions necessary to produce this 
evolution. Such models can display arbitrarily rapid expansion in the near future, leading to the 
destruction of all bound structures (a "little rip"). We determine observational constraints on these 
models and calculate the point at which the disintegration of bound structures occurs. For the same 
present-day value of w, a big rip with constant w disintegrates bound structures earlier than a little 
rip. 



I. INTRODUCTION 

Observations indicate that roughly 70% of the en- 
ergy density in the universe is in the form of an exotic, 
negative-pressure component, dubbed dark energy (HQ- 
(See Ref. [3] for a recent review.) If pus and pde are 
the density and pressure, respectively, of the dark en- 
ergy, then the dark energy can be characterized by the 
equation-of-statc parameter wde, denned by 



wde = Pde I Pde- 



(1) 



It was first noted by Caldwell Q that observational 
data do not rule out the possibility that wde < — 1- 
Such "phantom" dark energy models have several pecu- 
liar properties. The density of the dark energy increases 
with increasing scale factor, and both the scale factor 
and the phantom energy density can become infinite at 
a finite t, a condition known as the "big rip" [3-0] . It 
has even been suggested that the finite lifetime for the 
universe in these models may provide an explanation for 
the apparent coincidence between the current values of 
the matter density and the dark energy density Q. 

While w(a) < —1 as a extends into the future is a 
necessary condition for a future singularity, it is not suf- 
ficient. In particular, if w approaches —1 sufficiently 
rapidly, then it is possible to have a model in which pue 
increases with time, but in which there is no future sin- 
gularity. Conditions which produce such an evolution 
(specified in terms of pde as a function of pde) were 
explored in Refs. 0, 

In this paper, we examine such models in more de- 
tail. In particular, we will extend the parameter space 
discussed in Refs. [t| [l(| in both directions, showing 
that there are nonsingular models in which pe>e increases 
more rapidly than the nonsingular models discussed in 
those references, and, conversely, that there are singular 
models with pde increasing less rapidly than the singu- 
lar models discussed in Refs. [H, [lfj. Models without 
a future singularity in which pe>e increases with time 
will nonetheless eventually lead to a dissolution of bound 
structures at some point in the future, a process we have 
dubbed the "little rip." We discuss the time scales over 



which this process occurs. Finally, we consider observa- 
tional constraints on these models. 

In the next section, we examine the conditions neces- 
sary for a future singularity in models with w < — 1. In 
Sees. Ill and IV, specific little-rip models and disintegra- 
tion of bound systems are studied. Finally, in Sec. V, 
there is discussion. 



II. THE CONDITIONS FOR A FUTURE 
SINGULARITY 

We limit our discussion to a spatially flat universe, for 
which the Friedmann equation is 



P 
3' 



(2) 



where p is the total density, a is the scale factor, the 
dot will always denote a time derivative, and we take 
h = c = 8nG = 1 throughout. We will examine the fu- 
ture evolution of our universe from the point at which 
the pressure and density are dominated by the dark en- 
ergy, so we can assume p = pde and p = pde, and for 
simplicity we will drop the DE subscript. Then the dark 
energy density evolves as 



/> = -3( - ) (p + P ). 



(3) 



The simplest way to achieve w < — 1 is to take a scalar 
field Lagrangian with a negative kinetic term, and the 
conditions necessary for a future singularity in such mod- 
els have been explored in some detail jll] - [l4| . Here, how- 
ever, we explore the more general question of the condi- 
tions under which a dark energy density that increases 
with time can avoid a future singularity, and the conse- 
quences of such models. 

One can explore this question from a variety of starting 
points, by specifying, for example, the scale factor a as a 
function of the time t (an approach taken, for example, in 
Refs. [15l4l8| |). Alternately one can specify the pressure 
p as a function of the density p, as in Refs. [HG3. Note 
that this is equivalent to specifying the equation-of-statc 
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parameter toasa function of p, since w = p/p. Finally, 
one can specify the density p as a function of the scale 
factor a. Since we are interested specifically in nonsingu- 
lar models for which p increases with a, we shall adopt 
this last approach, but we will briefly examine the other 
two starting points. Of course, given any one of these 
three functions, the other two can be derived uniquely, 
but not always in a useful form. 

For example, suppose that we specify a(t). In order 
to avoid a big rip, it is sufficient that a(t) simply be a 
nonsingular function for all t. Writing 



J{t) 



(4) 



where f(t) is a nonsingular function, the density is given 
by equation © as p = 3 (a/a) 2 = 3/ 2 , and the condition 
that p be an increasing function of a is simply dp/ da = 
(6/d)// > 0, which is satisfied as long as 



/>0. 



(5) 



Thus, all little-rip models are described by an equation 
of the form with nonsingular / satisfying equation 

©• 

Now consider the approach of Refs. [3, [HI, who ex- 
pressed the pressure as a function of the density in the 
form 



V = -p- f(p), 



(6) 



where f(p) > ensures that the p increases with scale 
factor. In order to determine the existence of a future 
singularity, one can integrate equation ([3]) to obtain Q 



a = ao exp 



dp 

Wp)J' 



and equation ([2]) then gives [§} 



dp 



Zp-f(p)' 



(7) 



(8) 



The condition for a big-rip singularity is that the integral 
in equation (J8J) converges. Taking a power law for /(p), 
namely 



f(p) = A P a , 



(9) 



we see that a future singularity can be avoided for a < 
1/2 [UlltJ. We examine this boundary in more detail be- 
low, noting that one can have f{p) increase more rapidly 
than p 1 / 2 without a future singularity. 

Now consider the third possibility: specifying the den- 
sity p as an increasing function of scale factor a. We will 
seek upper and lower bounds on the growth rate of p(a) 
that can be used to determine whether or not a big-rip 
singularity is produced. Defining x = In a, we can rewrite 
equation ^ as 



(10) 




and the condition for avoiding a future big-rip singularity 
is 



1 



o \fpJx) 



dx — > oo. 



(11) 



The case p = —p — Ap 1 ! 2 from Refs. [TJ, [T(| corresponds 
to 



P_ 

Po 



3A 



ln(a/a ) + 1 



(12) 



where w < — 1 requires A > 0, and we take p = po and 
a = ao at a fixed time to. Expressing this density as a 
function of time rather than scale factor gives a much 
simpler expression: 



P_ _ p V3A(t~t ) 

Po 



(13) 



The equation-of-state parameter w corresponding to 
equation (TT2"j) can be derived from the relation 
(a/p)(dp/da) = -3(1 + 10): 



w 



-1 



1 



2 V a 1 ' A 

and the corresponding expansion law is 

2_ = e (2 v ^/3A)[ e (^ !A / 2 )< t - t o)-i] 
ao 



(14) 



(15) 



However, we can find p(d) for which p increases more 
rapidly with a, but for which equation (|11[) is still satis- 
fied. For example, writing p 1 ' 2 ~ (lna)(lnlna) as a — > oo 
satisfies equation (jlip . An example of such a p, with a 
free parameter B, is 

p_ = N f± B \ (i + H^ + B)) 2 (Hi + H^ + B))) 2 

Po W ) (l + Ml + B)) 2 (ln(l + In(l + B)))2 

(16) 



where the choice 



v \±b\= { ^ + B)2 

' ao' J (l + B) 2 (±) 2 



(17) 



leads to a real, nonnegative p and an analytic form for 
the behavior of ait): 



a ( ln(l + ln(l + B))« 

a 



/P0/3(t-tQ) 



(l + B)(l + ln(l + B)) ln(l + ln(l + B)) 



-V-B. 

This argument can be extended further. In general, if 
we denote ln^a;) = lnlnln .... ln(x), where the logarithm 
on the right-hand side is iterated j times, then any func- 
tion of the form 



p — (lna) 2 (ln 2 a) 2 (ln 3 a) 2 ...(ln m a) z 



(19) 
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satisfies equation (fTTj) as a —> oo and avoids a big-rip 
singularity. A density increasing as in equation (1191) leads 
to an expansion law of the form 



exp(exp(exp ...(exp(i))...)), 



(20) 



where there are m+1 exponentials. We have omitted the 
constants in equations ([T9|) and ([20|) for the sake of clar- 
ity. Equation ()20|) , while growing extraordinarily rapidly, 
is manifestly nonsingular. While an expansion law of this 
sort might seem absurd, it is probably less so than a big- 
rip expansion law, and in any case our goal is to try 
to determine the boundary between little-rip and big-rip 
evolution for p(a). In this spirit, consider the slowest 
growing power-law modification to equation (|19|) : 



p~ (lna) 2 (ln 2 a) 2 (ln 3 a) 2 ...(ln m a) 



2 + £ 



(21) 



where e > is a constant. No matter how small e is, 
and despite the fact that it modifies an extraordinarily 
slowly growing nested logarithm function, the growth law 
in equation (|21[) leads to a future big-rip singularity. 

Note that the bounds specified by equations (fll)j) and 
pTjl are not sharp; we can always find forms for p(a) that 
interpolate between these two behaviors and produce ei- 
ther a little rip or a big rip. However, as we take m to 
be arbitrarily large, nearly any function of interest will 
increase more rapidly than equation (|19j) or more slowly 
than equation ([2~Tjl . allowing us a practical, if not a rig- 
orously sharp, bound. This lack of a sharp bound is due 
to the fact that there is no bound on the fastest growing 
function a(t) which is nonsingular at finite t. 

If one is willing to place other restrictions on the form 
of p(a), then more stringent bounds apply. Barrow [l9| 
demonstrated that if p + 3p is a rational function of a 
and i, and a(t) is nonsingular at finite i, then a(t) can 
grow no more rapidly than the double exponential of a 
polynomial in t. Our equation (|19p violates this condition 
because of the logarithmic functions. 



III. CONSTRAINING LITTLE-RIP MODELS 

Here we shall examine in more detail the two specific 
little- rip models given by equations (|T2"j) and (IT51) , which 
we will call model 1 and model 2, respectively. Note that 
we do not make use of equations (fl~5f and (fl8|) here, as 
these are valid only when the matter density can be ne- 
glected in comparison to the dark energy density. Model 
1 is characterized by a single free parameter A, and the 
scale factor behaves asymptotically as a double exponen- 
tial in t, as in equation (|15[) : 



a[t) — > e 



(22) 



The parameter A is chosen to make a best fit to the latest 
supernova data from the Supernova Cosmology Project 
;20], and has the best-fit value A = 3.46 x 10~ 3 Gyr _1 , 
while a 95% C.L. fit can be found for the range —2.74 x 
lO^Gyr" 1 < A < 9.67 x Kr 3 Gyr -1 . 



Model 2 is characterized by the free parameter B and 
has a scale factor that behaves asymptotically as a triple 
exponential in t, as in equation (fT8|) : 



o(t) 



(23) 



The parameter B is chosen to make a best fit to [2(| 
as well, and it has the value B = 1.23. The confidence 
interval for B at the 95% C.L. is 1.12 < B < 1.34. In 
fitting both models, il mg — 0.274, Q Ko = 1 — fi mo , and 
Hq = 70.1 km s _1 Mpc -1 , which are consistent with the 
best-fit ranges for these values given by WM AP [2l[ . The 
resultant Hubble and residual ACDM (w = —1) plots of 
distance modulus p versus redshift z for both models are 
displayed in Fig. 1. 




FIG. 1: Top panel: Hubble plot of distance modulus fi ver- 
sus redshift z for the ACDM (w — —1) model (green) and 
models 1 (brown) and 2 (red). The lines are essentially indis- 
tinguishable. Bottom panel: The ACDM model is subtracted 
from models 1 (brown) and 2 (red). The ACDM model is, by 
definition, represented by the Ap = axis. As can be judged 
by the size of the error bars of the data, all are excellent fits 
to the supernovae data. (The color plots are in the online 
version of the paper.) 



Not surprisingly, the best-fit models closely resemble 
the ACDM model, which is known to be an excellent fit 
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to the data [22J. To see this more clearly, note that our 
models will resemble a cosmological constant at low red- 
shift as long as p(a) ~ constant for a ~ a . For model 1, 
this condition is satisfied when A/ y/po *C 1 in equation 
(p~2|) . while for model 2, we require B ~ 1.39 in equation 
(fl"6| . To see that i? should be close to this value, one 
should expand equation (|16p around a = a^. The zeroth- 
order term is poi and the coefficient for the first-order 
term is when B = 1.39. A comparison with our best-fit 
values indicates that these conditions are, indeed, satis- 
fied. Furthermore, in the limit where these conditions 
arc satisfied, these little-rip models closely resemble, at 
low redshift, big-rip models close to ACDM, i.e., mod- 
els with constant w < — 1 and |1 + w\ -C 1. To see 
this, recall that constant-w big-rip models have a density 
varying with scale factor as 

p = Po(a/a )-^ 1+w \ (24) 

For |1 + w\ <C 1 and a/ag not too far from 1, equation 
(f2~4"| behaves as 

p«p [l-3(l + u))ln(a/a )]. (25) 

Equation ([12]) reduces to equation flU} for A/^/po -C 1, 
with A/ y /p^= -(1+w). 



It is straightforward to estimate the corresponding 
tQ-($ for big-rip models with constant w to be [6j 

■~*(irf§r)- < 28 » 

and it is almost identical to t r i p , which is about one year 
later. 

Clearly, little-rip models can produce this disintegra- 
tion either earlier or later than big-rip models, depending 
on the exact parameters of each model. For example, by 
putting, w = — 1 — 1CP 3 in equation (J2SJ), we find a value 
of 11 Tyrs for i©-®, which is larger than that of models 1 
and 2 in Eqs. ([2l)l |2"7| . In this case, disintegration occurs 
earlier in the littlc-rip model than in the big-rip model. 

The five energy conditions (weak, null, dominant, null 
dominant, strong) (see, e.g., Ref. [24| ) are all violated 
by all little-rip and big-rip models. A simple way to 
see this is that if w < — 1, which occurs for any rip, a 
boost is allowed with (v/c) 2 > —w/c to an inertial frame 
with negative energy density. Having said that, if general 
relativity itself fails for length scales bigger than that of 
galaxies, we may not be constrained by the same energy 
conditions. 



V. DISCUSSION 



IV. DISINTEGRATION 

A feature of a big rip is that all bound-state systems 
disintegrate before the final singularity |5j . Here we show 
that little-rip models, despite not having a final singular- 
ity, also produce the disintegration of bound structures. 
As a first approximation, the disintegration time is when 
the dark energy density equals the mean density of the 
system. A more accurate method was presented in Q. 
We shall employ both methods to estimate the disintegra- 
tion of the Sun-Earth system. 1 For the little-rip models 1 
and 2, with the best-fit parameters derived in the previ- 
ous section, we find the time i©-® from the present time 
t until the Earth (©) - Sun (©) system is disintegrated 
to be: 

Model 1 : i©_© ~ 8 Tyrs (26) 

Model 2 : i Q _ ffi ~ 146 Gyrs. (27) 

Note that the disintegration time for model 2 is less than 
that of model 1, which is expected since p for model 2 
grows faster than p for model 1. 



1 When the Sun becomes a red giant in ~ 5 Gyrs, it will envelope 
Mercury and Venus, and (maybe) Earth |23| . Here, for the sake 
of making a point, we assume the Earth will continue to orbit 
the Sun until unbound by dark energy. 



In the big rip, the scale factor and density diverge in a 
singularity at a finite future time. In the ACDM model, 
there is no such divergence and no disintegration because 
the dark energy density remains constant. The little rip 
interpolates between these two cases; mathematically it 
can be represented as an infinite limit sequence which 
has the big rip and the ACDM model as its boundaries. 
Such models can be represented generically by a density 
varying with scale factor as in equation (|19|) . 

Physically, in the little rip, the scale factor and the 
density are never infinite at a finite time. Nevertheless, 
such models generically lead to structure disintegration 
at a finite time. For models consistent with current su- 
pernova observations, such disintegration can occur ei- 
ther earlier or later in a little-rip model than in a big- 
rip model, depending on the parameters chosen for the 
models. However, for a given present-day value of w, the 
big-rip model with constant w will necessarily lead to an 
earlier disintegration than the little-rip model with the 
same present-day value of w. This results from the fact 
that w increases monotonically in the little rip models, 
resulting in a smaller value for p at any given a than in 
the corresponding constant-w big-rip model, and there- 
fore, a lower expansion rate. Thus, supernova bounds on 
the epoch of disintegration for constant-u; big-rip mod- 
els also apply to little-rip models; one cannot simultane- 
ously satisfy supernova constraints and hasten the onset 
of disintegration to an arbitrarily early time simply by 
iterating exponentials in the expansion law. 

Furthermore, supernova data force both big-rip and 
little-rip models into a region of parameter space in which 



5 



both models resemble ACDM. In this limit, big-rip and 
little-rip models produce essentially the same expansion 
law up to the present, despite having very different future 
evolution. Thus, current data already make it essentially 
impossible to determine whether or not the universe will 
end in a future singularity. 

Finally, we remark that since the novel and specula- 
tive cyclic cosmology proposed in Ref. [25| requires only 
disintegration and not a singularity, such cyclicity would 
seem to be possible within a little-rip model instead of 
the big rip considered in [2^]. This is one potentially 
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